Chaos in Andreev Billiards 
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A new type of classical billiard — the Andreev billiard — is investigated using the tangent map tech- 
nique. Andreev billiards consist of a normal region surrounded by a superconducting region. In 
contrast with previously studied billiards, Andreev billiards are integrable in zero magnetic field, 
regardless of their shape. A magnetic field renders chaotic motion in a generically shaped billiard, 
which is demonstrated for the Bunimovich stadium by examination of both Poincare sections and 
Lyapunov exponents. The issue of the feasibility of certain experimental realizations is addressed. 
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In the development of the understanding of chaos, a 
prominent role has been played by the class of classical 
mechanical systems known as billiards Q. In such sys- 
tems, particles are confined by a step-like, single-particle 
potential to a region of space within which they propa- 
gate ballistically. Unless the shape of the billiard is highly 
regular (e.g., circular), in which case the system is inte- 
grable, the motion of a particle is chaotic. Unpredictabil- 
ity, the hallmark of chaotic motion, can be diagnosed 
qualitatively by the morphologies of Poincare sections, 
and more quantitatively by the corresponding Lyapunov 
exponents, the positivity of which signals the exponential 
sensitivity of trajectories to initial conditions. 

A common feature of all versions of billiards studied to 
date 0] is that reflection at boundaries is specular, i.e., 
only the component of the velocity normal to the bound- 
ary is inverted. We refer to such versions as conventional 
billiards (CBs; see Fig. |l|a, left). The purpose of this 
Letter is to explore the issue of classical chaos in a novel 
class of billiards, which have the property that scatter- 
ing at boundaries is retro-reflective, i.e., all components 
of the velocity are inverted, as is depicted in Fig. |ia 
(right). We refer to such billiards as Andreev billiards 
(ABs). Although we are unaware of examples of such a 
reflection mechanism in the realm of classical physics, a 
well-known example exists in condensed matter physics: 
Andreev reflection of electronic quasiparticles (having en- 
ergies in the superconducting gap A) from the normal-to- 
superconductor interface Thus, we envisage ABs as 
normal (N) domains surrounded by superconductor (S). 
It is adequate to regard the motion of electronic quasi- 
particles as semiclassical, provided that the billiard size 
is much larger than their typical de Broglie wavelength. 

The change from specular to Andreev reflection has 
a striking consequence in the context of chaos: whereas 
typical motion in a generically-shaped CB is chaotic, mo- 
tion in ABs is integrable, regardless of the shape of the 
billiard. This integrability becomes evident from the 



observation (cf. Fig. |l|a) that all motion occurs along 
chordal trajectories connecting only two points on the 
boundary: what system could be less ergodic? 

The presence of a magnetic field B substantially alters 
the situation (see Fig. |l|b), giving rise to a Lorentz force 
that curves the trajectories rit) of quasiparticles accord- 
ing to the equation of motion = {q/m)r x B. Now, An- 
dreev reflection inverts the quasiparticle charge q, mass 
m and velocity r, so that in the vicinity of the reflection 
point the acceleration of the outgoing hole is opposite to 
that of the incoming electron. Therefore, in a magnetic 
field the hole trajectory (dashed line in Fig. [l[b) no longer 
retraces the electron trajectory (full line) 13] , and vice 
versa, thus allowing the motion to explore the billiard. 
This raises the possibility of chaotic motion, a possibility 
that we explore in this paper. Our primary conclusion is 
that although in zero B ABs are integrable regardless of 
their shape, integrability is destroyed by the application 
of magnetic field for all but highly regular shapes. 

The study of ABs may provide an interesting link be- 
tween two rapidly developing fields: mesoscopic chaos 
Q and mesoscopic superconductivity Q . The prominent 
virtue of ABs, viz., that they are integrable in zero B 
regardless of shape and are rendered chaotic by nonzero 
B, makes them attractive from the experimental point 
of view. By comparison, the integrability of nanoscale 
CBs (e.g., in 2DEG heterostructures 1^) is immensely 
fragile, being readily destroyed by unintentional shape 
deformations or surface roughness. Therefore, to obtain 
nontrivial chaos in CBs, i.e., chaos caused by intentional 
choice of shape, one must use state-of-the-art nanofabri- 
cation technology. In contrast, an AB prepared without 
any special attention to shape will be integrable, nonin- 
tegrability of varying degree being achieved by adjusting 
an external parameter, viz., B 

To explore qualitatively the implications of a magnetic 
field for the integrability of ABs, we focus on a planar 
two-dimensional AB in a magnetic field perpendicular to 
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the plane. We take the dynamics to be classical cyclotron 
motion with radius Rc (= mv/qB) of a particle inside 
the billiard, supplemented by Andreev reflection at the 
boundary. The case of CBs in a magnetic field has been 
studied extensively and the tangent map PJlO| 

has proven to be a convenient approach. A tangent map 
is a variant of a Poincare map, in which the state of the 
system is monitored only at collisions with the boundary, 
the remainder of the motion being obtained by a simple 
geometry. Thus, the problem reduces to one of following 
the sequence of reflection points generated by the particle 
as it explores the billiard. 

A priori, our phase space is four-dimensional: two 
components of the position in the plane and two conju- 
gate momenta. Energy conservation constrains the mag- 
nitude of the momentum, leaving one freedom, which (as 
we are using the tangent map) we take to be the angle 
a between the velocity v and the tangent to the bound- 
ary (see Fig. In addition, the fact that reflections 
take place on the boundary leaves one further freedom, 
which we take to be the arclength s along the boundary. 
Then the (continuous time) dynamics is replaced by the 
(discrete) map: {s,a} {s'(s, a), q;'(s, a)}, embodying 
cyclotron motion followed by Andreev reflection. For the 
sake of convenience, we monitor only reflections of quasi- 
particles of the same (say, electron) type. (Thus, e.g., for 
the sequence of reflections ^ 1 ^ 2 shown in Fig. |^a, 
only reflections and 2 are used to construct Poincare 
sections.) 

First, consider the case of an AB of arbitrary shape 
at i? = 0. In this case, all trajectories, such as that 
depicted in the right billiard of Fig. |^a, are trivially pe- 
riodic, and the Poincare section for a given trajectory 
reduces to a single point, completely determined by the 
initial conditions. We analyze the case of B for 
the example of an AB in the shape of a Bunimovich 
stadium |ll[], as shown in the top row of Fig. ||. Fig- 
ure ^ shows Poincare sections (bottom row) for a selec- 
tion of initial conditions, along with typical trajectories 
(top row). For convenience, we introduce the dimension- 
less magnetic field f3 = R/Rc oc B and the tangential mo- 
mentum p = cos a. For the case of a weak field {f3 — 0.02, 
left column), {s,p} space is apparently foliated by well- 
defined curves, each curve corresponding to a particular 
initial condition. Although it appears that the motion 
is integrable, when viewed at a finer scale one sees the 
breakdown of foliation, as shown in the inset. Thus, 
the motion is in fact weakly chaotic, as we have also 
confirmed by examining the corresponding Lyapunov ex- 
ponent. In intermediate fields (for which the cyclotron 
radius is comparable to the billiard size), the Poincare 
sections (except those in a central region) appear to fill 
a two-dimensional area of the {s,p} plane with discon- 
nected points ("dust"), as shown for the case /? = 0.33 
in the middle column. Such behavior is commonly taken 
as an indication of chaos and thus our Poincare 

sections suggest that indeed ABs are rendered chaotic 
by the application of magnetic field. Notice the scar-like 



feature running across the {s,p} plane: it is a remnant 
of the foliation that dominates in weaker fields. A sin- 
gle, additional, quasi-one-dimensional Poincare section, 
corresponding to an initial condition deliberately chosen 
in the scar, is also shown. In strong fields, when the 
cyclotron radius is much smaller than the billiard size, 
particles move along skipping trajectories. On the scale 
of a typical skip, there is little distinction between motion 
over straight and semi-circular segments of the boundary 
and, therefore, the motion is less sensitive to the bil- 
liard shape. As a result, the Poincare section exhibits 
a (partial) re-entrance of integrability, i.e., the "dust" 
that arises in intermediate fields is reorganized into some 
structures, as is seen in the bottom row (for /3 = 10). 
This structure bears a certain resemblance to that found 
in weak fields, thus indica ting a trend towards less chaotic 
behavior. As with CBs |10(] , chaos is most pronounced 
for intermediate fields, becoming less pronounced in both 
the weaker and stronger field regimes. 

To provide quantitative support for the suggestion of 
chaos inferred from the inspection of Poincare sections, 
we now turn to the computation of Lyapunov exponents, 
which characterize the rate of exponential divergence 
of trajectories having initial conditions nearby in phase 
space. This is accomplished by investigating the stabil- 
ity of the tangent map via the adaptation to ABs of the 
method of Rcfs. 10|. Consider the situation depicted in 
Fig. |l|c. From the kinematics of circular motion we have: 
Vq — LJc X 'I'l = vq — X ro, where ujc = qB/m and 
rQ i are the radius- vectors of the reflection points; other 
notations are defined in Fig. The tangent map is de- 
rived by varying this equation with respect to s and p, 
and relating the deviations of two nearby trajectories (5q 
(= {ds,Sp}) before ((5qo) and after ((5qi) refiection from 
point 1: Sqi = 7i_o • Sqf). After some straightforward 
algebra, we find that Ti^ is given by 
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where po,i are the radii of curvature of the billiard bound- 
ary at the reflection points and 1, (j{(f>) = sin0, and x is 
defined in Fig. |l|c |l2j. After N bounces, the separation 

Jqat is determined by the matrix TJv,o = IljLi'^j-i- 
The largest Lyapunov exponent associated with a given 
trajectory is calculated as A = lim^v^oo ^n, where 



TrrAr,o/2 + V(Tr Vo/2)' - 1 ■ (1) 



We have calculated the Lyapunov exponents for a wide 
selection of initial conditions {so,Po} and values of B. A 
typical sequence Xn is shown in Fig. ^. The convergence 
of Aat to a non-zero value as ^ oo provides evidence 
for the exponential divergence of nearby trajectories, i.e., 
chaos. 

We now turn to the issue of possible experimental real- 
izations of ABs. In one possible scheme, an AB is formed 
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by surrounding a 2DEG with a superconducting contact 
||13| . The chaotic nature of the motion can be diagnosed 
either by passing normal current through the structure 
and measuring the conductance, as with nanoscale CBs 
, or by studying, e.g., via STM, correlations in real- 
and energy-space, which provide signatures of classical 
chaos at the quantum level jist . The primary demand 
on the experimental realization of all billiards, includ- 
ing ABs, is that the motion of the electrons inside the 
billiard be ballistic. Our analysis of Poincare sections 
and corresponding Lyapunov exponents shows that when 
L ^ R, chaos is established after a few bounces, 
and thus it will not be masked by impurity scattering 
provided that L,i? ^ i^, where 1^ is the elastic mean 
free path of the N region. On the other hand, B should 
not exceed the (lower) critical field of the superconductor 
Be and, hence, R^ > i?""" = pp/eBc, where pp refers to 
the N region. Thus, it is sufficient to have £c ^ i?™™. 
Taking parameters for the Nb/InAs structure studied re- 
cently (density of electrons Uc — 9 x 10^^ cm^^, Ref. [ p^ ; 
Bf. Ki 2000 G) we obtain £e > 0.6 ^m, which is ac- 
cessible via current nanofabrication technologies. (E.g., 
£e ~ 3.5 /im in Ref. 10.) A drawback of the scheme de- 
scribed above is that the superconductor and the 2DEG 
are metallurgically distinct, and thus the probability for 
normal scattering at the interface is nonzero, at the ex- 
pense of the Andreev reffection, which results in billiard 
having a mixed AB/CB character. This drawback can 
be eliminated by employing the proximity effect. The 
AB is formed in a region of a superconductor where the 
superconductivity has been suppressed, either due to the 
vicinity of a normal metal island (see Ref. ) , or (with 
a type I superconductor) by the application of a mag- 
netic field, which creates domains of N phase. As such 
a scheme eliminates metallurgical boundaries between N 
and S, reflection is of purely the Andreev type. 

An interesting direction of further research would be 
to explore the quantum mechanics of ABs. At least three 
directions are immediately apparent: (i) spectral geome- 
try and the Weyl-Kac problem; (ii) energy-level statis- 
tics, random matrix approaches and universality; and 
(iii) spatial structure of quasiparticle wavefunctions. 
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FIG. 1. (a) Typical trajectories for specular (SR) and 
Andreev (AR) reflection at B = 0. The CB (left) is formed by 
a single-particle potential, the AB (right) by a pair-potential, 
(b) AR and SR in a magnetic fleld. For AR, each cyclotron 
orbit is necessarily tangential to the previous one, in contrast 
with SR. (c) Geometry of the tangent map. A particle starts 
from with velocity vo, follows the cyclotron trajectory across 
the billiard, arriving at 1 with velocity Vq. Owing to the 
nature of AR, the velocity vi after the reflection at 1 is — Vq. 
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FIG. 2. Top row: a typical trajectory for a Buni- 
movich-stadium-shaped AB {L = R) for 3 values of the mag- 
netic field: (a) /? = 0.02; (b) (3 = 0.33; (c) /? = 10. Bottom 
row: Poincare sections in these fields constructed by follow- 
ing the first 1000 bounces for the trajectories starting with 
ao = 10°, 20°, ... , 170° (cf. Fig. |l|c) from random points on 
the perimeter of the billiard. Thin vertical lines on Poincare 
sections separate regions corresponding to straight (wider) 
and semi-circular (narrower) segments of the billiard bound- 
ary. In the Poincare sections for the weak field, flat segments 
result from almost-chordal motion across a single semicir- 
cle. {Only flat regions would occur for a circular billiard.) 
Similarly, curved regions result from trajectories connecting 
any two of the four distinct segments of the boundary. In- 
set: segment (indicated by arrow) of the foliation, magnifled 
(a;x ~ 10; yx ~ 50). Ticks on the stadium boundaries mark 
the points s — 0; flUed circles indicate the start of trajectories 
So- We choose units in which the billiard perimeter is unity. 

FIG. 3. Lyapunov functions Ajv [log-scale; see Eq. (^)] vs 
number of bounces N for: (b) (3 = 0.33; (c) /3 = 10. The 
associated trajectories are shown at the top of Fig. |^. For 
13 = 0.02 (not shown) we find A ^ 0.002. 



